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We present a formalism able to predict the transformation of light beams passing through biaxial 
crystals. We use this formalism to show both theoretically and experimentally the transition from 
double refraction to conical refraction, which is found when light propagates along one of the optic 
axes of a biaxial crystal. Additionally, we demonstrate that the theory is applicable both to non- 
cylindrically symmetric and non-homogeneously polarized beams by predicting the transformation 
of input beams passing through a cascade of biaxial crystals. 
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I. INTRODUCTION 

Light propagation in optically anisotropic media is an 
interesting object of study due to the number of scien¬ 
tific and technological applications of uniaxial and biax¬ 
ial crystals. Although the effects of anisotropic media on 
light beams is known since centuries ago, even nowadays 
this topic is an object of exhaustive study. Due to the fact 
that the mathematical description of uniaxial crystals is 
simpler than for the case of biaxial crystals, light propa¬ 
gation along the former is well understood. However, to 
our knowledge there is no basic set of equations offering 
the possibility of predicting the behavior of light after 
passing through biaxial crystals along any beam propa¬ 
gation direction. The most related works on that topic 
have been performed by Dreger [l| and Gamier [ij . Anal¬ 
ogously to uniaxial crystals, in biaxial crystals light prop¬ 
agating out of one of the optic axes suffers from double 
refraction. In contrast, if the beam propagates along one 
of the optic axes it suffers from conical refraction (CR), 
being this a phenomenon of renewed interest during the 
last years. In this case, the transformation of cylindri- 
cally symmetric and homogeneously polarized beams has 
been well described with diffractive optics by Belsky and 
Khapalyuk, and Berry gi- Under conditions of CR, 
an input Gaussian beam with waist radius wq, is trans¬ 
formed at the focal plane of the system into a pair of 
concentric bright rings split by a dark (Poggenforff) ring 
of geometrical radius Rq provided that po = Ro /'^o ^ 1 
(see Fig. Hd)). No matter what the state of polariza¬ 
tion (SOP) of the input beam is, at any point of the 
rings the SOP is linear with the azimuth rotating along 
the rings so that every two diametrically opposite points 
have orthogonal polarizations [5|. Out of the focal plane, 
the beam evolves as an op tical bottle with two on-axis 
maxima at positions 2 : = @ (where zr is the 

Rayleigh length of the Gaussian input beam). 


* Corresponding author: alejandro.turpin@uab.cat 


Up to know it has been shown both experimentally 
and theoretically CR for homogeneously polarized Gaus¬ 
sian input beams diQ, tcm-hat input beams (H and 
Laguerre-Gauss input beams [9|,[l3 through a single crys¬ 
tal and also the transformation of Gaussian input beams 
in cascaded crystals 0, liMl. With a wave-vector 
based formalism, the CR transformation of elliptical in- 
put beams and an axicon input beam has been analyzed 
[l5| as well. Berry and Jeffrey also showed the transition 
from double refraction to CR for cylindrically symmet¬ 
ric beams with a modified treatment of CR that allows 
studying dichroic crystals M- 

Our aim here is to present a theoretical formalism ca¬ 
pable to predict the transformation of light beams in 
biaxial crystals independently on the beam shape, its 
SOP and its propagation direction. We will provide 
the formulas for the electric field behind the crystal and 
show the evolution from double refraction to CR and the 
transformation of non-cylindrically symmetric and non- 
homogeneously polarized beams both in single and mul¬ 
tiple biaxial crystals. As a proof of the usefulness of the 
theory, we will present the experimental cases of a Gaus¬ 
sian and an elliptical input beam. 


II. THEORETICAL FORMALISM 

In what follows, we will consider normalized coordi¬ 
nates to Wo and zr, i.e. x xjwo^ y -T y/wo and 
z ^ z/ Zr. In the parabolic approximation, after passing 
through a medium or optical element, a light beam can 
be described by means of its displacement vector D as 
a superposition of plane waves R = which are 

generated from a unitary transformation provided by the 
optical element U(R) applied over the Fourier transform 
vector of the input light beam A{R). In other words, 
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D =JJ dK:,dKye^'^'^U{fl)A{K), ( 1 ) 

— 00 
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where A{k) is the Fourier transform 

- Axi^lAjCx Ayi^lAjCy ^ ( 2 ) 

oo 

Ax{k) = Jj Ex{f)e~'‘^'^dxdy, (3) 

— OO 
OO 

^y{«) = Jj Ey{f)e-^^-^dxdy, (4) 

— OO 

of the input beam with transverse electric field 

E{r ) = Ex^x + EyCy^ (5) 


where 03 and <Ji are the Pauli matrices and I is 2 x 2 
the identity matrix. It is straightforward to demon¬ 
strate that for a given vector v = vn with |n| = 1 , 
the evaluation of gives Icos(u) + i{n • asin(v))^ 

where a = (cfi, cf 2 , 03 ) is the vector of Pauli matrices. 
By recalling the latter identity and after evaluation of 
by using Eq. (| 6 ]), the unitary transforma¬ 
tion provided by the crystal can be obtained: 


where and are the unitary vectors in Cartesian co¬ 
ordinates. 

For a low birefringent biaxial crystal, it has been shown 
that the unitary transformation provided by the medium 
is tf {k) = with [ 3 , 0 , [l 3 | 

T{fi) = +poR- ( 6 ) 

_I 


xf , .sin(poK) ^ / 

cos{/9oK)I + ^-K-{a-3,cri) , 

(7) 

where k = ^+ Ky. By combining Eq. ([T]) with Eq. ([Tj), 

there can be obtained two main integrals that describe 
the beam evolution behind the biaxial crystal: 


U{k) = e“*2'- 


Bo^a{r,po) 

Bi,a{r,po) 




2n'^ Sm{pQK,) Aa{K)dK,xdK,y, 


— 00 
00 
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^ II ^ (cos(po«;) 


+ i — sin (pok) ) Aa{i^)dKxdKy, 

K, ' 


(8) 

(9) 


where a = {x,y) and n is the mean refractive index of 
the biaxial crystal. The expressions for the field D in 
terms of Eqs. (jS]) and ^ are 

Dx = Bo^y{f, po) + Bi^x{e Po), (10) 

Dy = Bo^x{ePo) E Bi^y{f,-po). ( 11 ) 

Note that Eqs. ©-dm) do not require the input 
beam to be cylindrically symmetric nor to be homoge¬ 
neously polarized. Additionally, since an unpolarized 
beam can be expressed as the incoherent addition of mul¬ 
tiple beams, these equations can be used to predict the 
transformation of any input beam, no matter its SOP 
or its shape, as long as its Fourier transform can be ob¬ 
tained. For input beams with non-homogeneous polar¬ 
izations, it must be taken into account that the beam 
and the polarization can be always decomposed in the 
(x, y) basis, so that the theoretical formalism presented 
above can be always used. The above formalism reduces 
to the well known Belsky-Khapalyuk-Berry equations 
for CR when cylindrically symmetric and homogeneously 
polarized beams are considered, as demonstrated in Ap¬ 
pendix 

In what follows we will demonstrate the flexibility of 
the presented formalism to show the transformation of 


light beams in a variety of situations. 

III. SINGLE CRYSTAL CONFIGURATION 

In this section we discuss the transformation of in¬ 
put beams that propagate within a biaxial crystal non¬ 
parallel to one of the optic axes, i.e. under conditions 
of double refraction. We consider homogeneously left 
handed circularly polarized input beams, i.e with nor¬ 
malized electric field eo = (l,i)/\/ 2 , and we look at the 
transverse patterns at the focal plane {Z = 0). Propaga¬ 
tion out of the optic axis can be described by the prod¬ 
uct of the input field amplitude by a factor exp(—i^-fA), 
where 6 = (Sx^Sy) and r± = (x^y) 0 . 

We consider the well known case of an input beam with 
Gaussian transverse profile and an elliptical input beam. 
Their electric field amplitudes are described by 

Eg{x, y) = (12) 

EB{x,y) = (13) 

where a and b are constants and Sx and Sy give the angu¬ 
lar separation of the input beam’s propagation direction 
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with respect to the optic axis of the crystal. 

The experiments have been performed by employing 
a KGd(W 02)4 (conicity a = 16.9 mrad) biaxial crystal 
with a length of / = 28 mm, yielding a CR ring radius of 
Ro = la = 475 /rm. In FigslT](bl)-(b5) we used a Gaus¬ 
sian input beam focused by a spherical lens with 100 mm 
of focal length. The elliptical input beam was obtained 
by focusing the same Gaussian beam with a cylindrical 
lens with 100 mm of focal length. The biaxial crystal 
was mounted on an angular micrometric positioner that 
allowed to change the 0 and 0 angles in spherical coor¬ 
dinates and to observe the transition from double refrac¬ 
tion to GR as the optic axis and the beam propagation 
direction approached each other. 

Fig. [ 1 ] shows transition from double refraction to GR 
for the Gaussian input beam (first and second rows) and 
the elliptical input beam (third and fourth rows) both 
experimentally (second and fourth rows) and numerically 
calculated by using Eqs. (|3l)- (pT]) (first and third rows). 
For simplicity, we consider only angular displacement of 
the input beam in the vertical direction, i.e. = 0 . 
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FIG. 1. Transverse intensity patterns showing the transi¬ 
tion from double refraction to conical refraction for a Gaus¬ 
sian input beam (first and second rows) and an elliptical in¬ 
put beam (third and fourth rows) both experimentally (sec¬ 
ond and fourth rows) and numerically calculated by using 
Eqs. @-([13]) (first and third rows). 

The transformation of a Gaussian input beam propa¬ 
gating parallel to one of the optic axes of a biaxial crys¬ 
tal as described by Eq. (parallel propagation im¬ 
plies Sx = = 0 ) is the most commonly studied case 

in GR. When Sx,y 7 ^ 0 double refraction instead of GR 
is found. Double refraction in uniaxial crystals is as¬ 
sociated with the splitting of the input beam into two 
beams with identical transverse pattern and orthogonal 
polarizations, corresponding to the ordinary and the ex¬ 
traordinary polarizations. However, in biaxial crystals, a 
Gaussian input beam propagating out of the optic axes 
splits into two azimuthal sectors placed at diametrically 


opposite positions of the otherwise expected GR ring, 
provided that the angular propagation deviation with re¬ 
spect to the optic axis is small. Only when 6x,y ^ 1 the 
output split beams preserve the input beam’s pattern, 
as in uniaxial crystals. As the beam propagation direc¬ 
tion approaches to the optic axis, the split beams occupy 
a larger azimuthal angle and eventually both interfere, 
see Eigs. [T]^a3)-(a5). The fact that the two split beams 
interfere implies that both beams possess regions of non- 
orthogonal polarizations. Eor parallel propagation with 
respect to the optic axes, the interference between both 
split beams is maximum and the two bright rings with 
Poggendorff splitting are formed. 

In case of an elliptical input beam, there is a competi¬ 
tion between the ellipticity of the shape of the beam and 
the refraction induced by the crystal. We consider an 
elliptical beam described by Eq. ([T3|) with a = 1 , 5 = 0 . 1 , 
i.e. with Wx < Wy. Since the misalignment with re¬ 
spect to the optic axis is only along the y direction, the 
two split beams are expected also to appear in that di¬ 
rection. Due to the non-symmetrical nature of the el¬ 
liptical beam and the double refraction provided by the 
biaxial crystal, which induce opposite effects, the two re¬ 
fracted beams for beam propagation out of the optic axis 
are wider than for the Gaussian input beam case, see 
Eig.[T](cl). As the misalignment of the input beam is re¬ 
duced, the refracted beams expand along the azimuthal 
direction and at some point both interfere, as in the pre¬ 
vious case. Eor parallel propagation to one of the optic 
axis, the pattern is formed by two lobes each of which 
with Poggendorff splitting, see Eig.[T]^c5). The two lobes 
are slightly connected between each other but we have 
checked that the connection points tend to disappear as 
the ratio of the axes of the ellipse increases. We have 
additionally checked that there is a continuous evolution 
of the double-concentric ring structure from Eig. [T]^a5) 
into the double-lobe pattern from Eig. HJcS) as the ratio 
of the axes of the ellipse increases. Eor a deeper study of 
GR from elliptical beams, see Ref. Hi- 


IV. CASCADE OF CRYSTALS 

Light propagation through a cascade of crystals can be 
predicted by using Eqs. ©-(HU recursively. In this case, 
the output electric field after the first crystal Ecru be¬ 
comes the input electric field impinging the second biax¬ 
ial crystal and this process can be repeated until passing 
through the whole cascade. Since the beams obtained 
after the first crystal can be non-cylindrically symmetric 
and non-homogeneously polarized, the results reported in 
what follows prove also the usefulness of our formulation 
for such input beams. Eig. [ 2 ] presents the numerical simu¬ 
lations and experimental results obtained for a cascade of 
two biaxial crystal with common optic axis being paral¬ 
lel to the optical axis of the system for an elliptical input 
beam. The passage of an elliptical input beam through 
the two crystal cascade shows a similar behavior: each of 
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the two lobes with Poggendorff splitting obtained after 
the first crystal split into two lobes, also with Poggen¬ 
dorff splitting, after passing through the second biaxial 
crystal. The experimental observations from Fig. EJb) 
were performed by using two KGd(W 02)4 biaxial crys¬ 
tals with lengths of h = 28 mm and I 2 = 10 mm. The 
Gaussian and elliptical input beams were identical to the 
ones used in the previous section. As it can be appreci¬ 
ated, both theory and experiment agree well. 


(a) (b) 

i I I i I I i 


FIG. 2. Transverse intensity pattern at the focal plane for 
a cascade of two biaxial crystals with aligned optic axes ob¬ 
tained from an elliptical input beam, (a) Numerical simula¬ 
tions obtained by using Eqs. (|8|- ([T3l) . (b) Experimental mea¬ 
surements obtained by using two KGd(W 02)4 biaxial crystals 
with lengths of li — 28 mm and I 2 = 10 mm. 


a wide variety of beams, including non-cylindrically sym¬ 
metric and non-homogeneously polarized input beams, 
propagating along one of the optic axes of a biaxial crys¬ 
tal. Additionally, we have shown the transition from dou¬ 
ble to conical refraction for non-cylindrically symmetric 
beams, such as an elliptical beam. Finally, we have re¬ 
ported the case of a cascade of multiple biaxial crystals. 
Our numerical predictions accurately reproduce previous 
experiments on light transformation in biaxial crystals 

diEsiii. 

We expect the reported results to be applied in all ar¬ 
eas of optics requiring the use of biaxial crystals and, in 
particular, in optical trapping [^, [T^-0L space opti¬ 
cal communications 0, Dolarimet rv 1^ . , lasing and 

non-linear optical phenomena l25]. laser processing 
(M, the generation and manipulation of vector beams 
U and beam shaping . 
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V. CONCLUSIONS 

In summary, we have presented a theoretical formalism 
that describes the beam evolution of light propagating 
through biaxial crystals. The method can be used with 
complete generality, i.e. it does not require symmetries 
in the beam shape nor in its state of polarization, as well 
as can be used to predict the behavior of light with any 
propagation direction. We have demonstrated the flexi¬ 
bility of the reported formalism for the transformation of 


Appendix A: Cylindrically symmetric solution 

In this Appendix we will demonstrate that the general 
formalism presented in Section [III can be used to derive 
the fundamental Belsky-Khapalyuk-Berry equations of 
conical refraction for input beams possessing both cylin¬ 
drical symmetry and homogeneous polarization. By us¬ 
ing Eq. ([7j), Eq. ([!]) can be written in Cartesian coordi¬ 
nates as follows: 


00 



—00 


^ .sin(poK) ^ ' 

COS(/)oK)I + l - K ■ (0-3, CTi) 


A(jV)dKxdKy. 


(Al) 


The electric field described by Eq. (EH) can be rewrit¬ 
ten in terms of two fundamental equations as shown by 

Eqs. (IHI-idID. 


If we change the coordinate system from Cartesian to 
cylindrical, we must carry out the following replacements: 


f -H- p{cos((fi), sin((^)), 

^ ^iKpcos((f)^-cp) 


(A2) 

(A3) 

(A4) 


// 


dKxdK,y ^ 


— 00 



(AS) 


where we consider the radial coordinate to be normalized 
to the beam waist radius (wq) of the input beam through 
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p = rlwo. In this case, Eq. ()A1D reads as follows: 


D = 


1 



^ikp cos{4)^-<p) Z‘^ 


cos(po/^) 


1 0 
0 1 


+ isin(po/^) 


cos(0^) sin(0;,) 

sin(0^) -cos((/)^) 


A{i^)hzdhvd(j)^ 


(A6) 


Now, we simplify the system by considering a uni¬ 
formly polarized and cylindrically symmetric input beam 
E = E(r) eo with cylindrically symmetric Fourier trans¬ 
form too, i.e. = Ay{R) = a(/^), where 


the following expressions 

2tv 

0 

27r 

J ^iK,pcos(cf)^-ip) 0 os( 0 /^)(i/T: = cos{(p)Ji{kp) 
0 


27T 


/ 


^iKpcos((f)^-(p) 


sm{(j)^)dhv = 0 


(A8) 

(A9) 

(AlO) 


CX) 

a{K) = J KE{r)Jo{rK)dr 


(A7) 


and obtain a ID integral for the D field: 

D ^ { ^0 + ^1 cos((^) Bi sin((/9) 

Bi sin((p) Bo — Bi cos{ip) 


eo, 


(All) 


where Bq = Bo{p, Z, po) and Bi = 5i(p, Z, po) are the 
Belsky-Khapalyuk-Berry integrals that describe the evo¬ 
lution of the CR beam both in the radial and axial di¬ 
rections: 


being Jn{x) the n^^-order Bessel function of the first type. 
With this assumption, we can integrate over by using 


oo 

1 r . ^2 ^2 

Bo{r,Z,po) = -— / a(K}e~'‘ cos{Kpo)Jo{Kr)^l^) 
(271-) J 


Bo{r,Z,po) = —— / a{K)e *'‘ 2 n sin(Kpo)Ji(7cr)((M3) 
(27r) J 
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